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Abstract 

The Bogoliubov-de Gennes equations are used for a number of theoretical works 
on the trapped Bose-Einstein condensates. These equations are known to give the 
energies of the quasi-particles when all the eigenvalues are real. We consider the 
case in which these equations have complex eigenvalues. We give the complete set 
including those modes whose eigenvalues are complex. The quantum fields which 
represent neutral atoms are expanded in terms of the complete set. It is shown that 
the state space is an indefinite metric one and that the free Hamiltonian is not 
diagonalizable in the conventional bosonic representation. We introduce a criterion 
to select quantum states describing the metastablity of the condensate, called the 
physical state conditions. In order to study the instability, we formulate the linear 
response of the density against the time-dependent external perturbation within 
the regime of Kubo's linear response theory. Some states, satisfying all the physical 
state conditions, give the blow-up and damping behavior of the density distributions 
corresponding to the complex eigenmodes. It is qualitatively consistent with the 
result of the recent analyses using the time-dependent Gross-Pitaevskii equation. 
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1 Introduction 



In 1995, the Bose-Einstein condensates (BECs) of neutral atoms in the mag- 
netic trap were realized [1,2,3]. After several years, some kinds of quantized 
vortices in condensates have been observed, e.g., singly quantized vortices and 
a vortex lattice [4,5,6]. Furthermore, the doubly quantized vortices have been 
created under the use of phase imprinting technique [7]. It has also been ob- 
served that a doubly quantized vortex decays into two singly quantized vortices 
[8]. 

Excitations around the condensates, with or without vortices, are represented 
by the solutions of the Bogoliubov-de Gennes (BdG) equations [9,10,11], which 
follow from linearizing the time- dependent Gross-Pitaevskii (TDGP) equation 
[12] and form a set of two-component eigenf unctions. If all the engenvalues 
are real, the corresponding free Hamiltonian is diagonalized in the bosonic 
representation and the quasi-particle picture is obtained. 

Recently, it is found that the BdG equations have complex eigenvalues in some 
cases, e.g., the case where the condensate flows in an optical lattice [13], or 
where the condensate has a highly quantized vortex [14,15,16,17,18] or gap 
solitons [19], or the case of the multi-component BECs [20,21]. It is also men- 
tioned that a subject of the complex eigenvalues is taken up in the context 
of diagonalization of the Hamiltonian of the quadratic form of creation and 
annihilation operators [22] and in the studies of the Gross-Pitaevskii (GP) 
energy functional [23]. These eigenfunctions with complex eigenvalues cause 
the "dynamical instability", which is a type of instability of the condensates 
and should be distinguished from the "Landau instability" caused by nega- 
tive energies of the quasi-particle. The former is associated with the decay of 
the initial configuration of the condensate, while the latter corresponds to the 
absence of thermodynamic equilibrium [12]. The dynamical instability of the 
condensates in optical lattices was observed in the experiment, and the values 
of the quasi-momentum at which the experimental loss rate has its peak are 
in good agreement with those at which the imaginary part of the complex 
eigenvalue of the BdG equations has its peak [24]. The time scale to form 
the ferromagnetic domains in the spinor BECs, observed experimentally by 
rapidly quenching the system to conditions in which the ferromagnetic phase 
is energetically favored, is roughly described by the imaginary part of the com- 
plex eigenvalue of the BdG equations [25]. In addition, when the condensate 
has a highly quantized vortex, the emergence of the complex eigenvalues is 
interpreted as the sign of dynamical instability caused by the external pertur- 
bations [17,18]. However, no corresponding experiment in which the conden- 
sate is perturbed by an external force has been reported, although the decay 
of a doubly quantized vortex into two singly quantized vortices, without per- 
turbation, is observed experimentally [8]. The relation between the complex 
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eigenvalues and the instability of highly quantized vortices is not elucidated 
fully and is still under study [26,27]. 

So far, the interpretation of complex eigenvalues in the BdG equations has 
been based on the TDGP equation, but not on the quantum field theory 
(QFT). It has not been made clear how the complex eigenmodes should be 
interpreted in QFT, and how to describe the dynamical instability against ex- 
ternal time- dependent perturbations within the formulation of QFT. In gen- 
eral, it is possible in terms of equilibrium QFT to give dynamical description 
in the linear regime, using Kubo's linear response theory (LRT) [28,29,30]. 
But Kubo's LRT is based on the canonical commutation relations (CCRs) 
and representations of the states. In the case where complex eigenvalues arise, 
it is not easy to keep the CCRs and to discuss the representations of states due 
to the indefinite metric which is characteristic of the BdG equations. It is well 
known in QFT that when the indefinite metric is involved, careful treatments 
on the complete set are needed (for example, the ghosts in gauge fields [31]). 
But this point has not been discussed in the framework of QFT with complex 
eigenvalues of the BdG equations and the meaning of the complex eigenmodes 
in Kubo's LRT has not been clear. 

In this paper, we construct QFT and formulate Kubo's LRT in the case where 
complex eigenvalues arise in the BdG equations. The blow-up or damping 
behavior of the density distributions, corresponding to the complex modes 
against the time-dependent external perturbation, is shown in the linear regime. 
It is consistent with the result from the recent analyses based on the TDGP 
equation [17,18]. Furthermore, we discuss that the effect beyond the GP the- 
ory is given if we consider the higher order terms of quantum fluctuations 
using the formalism established in this paper. 

This paper is organized as follows. In Sec. 2, the model action and Hamilto- 
nian are given. In Sec. 3, the BdG equations are introduced, and when all the 
eigenvalues are real, the mathematical properties of eigenfunctions including 
the zero-mode are summarized. In Sec. 4, for the sake of convenience in the fol- 
lowing discussions, the doublet notations are introduced. In Sec. 5, we consider 
the case in which the complex eigenvalues arise in the BdG equations. Mathe- 
matical properties of the eigenfunctions are summarized and the complete set 
of functions is given. In Sec. 6, the quantized fields are expanded using the 
complete set of functions given in Sec. 5. We check that the CCRs are kept in 
this expansion. Next, we study the properties of the Hamiltonian. Particularly, 
we show that the free Hamiltonian can not be diagonalized in the usual bosonic 
representation. Then we construct eigenstates of the free Hamiltonian, which 
are not normalizable. We discuss the complete set of states, corresponding to 
the complex eigenspace. In Sec. 7, we introduce the conditions for physical 
states, which are not only matched with the result of the experiment in which 
highly quantized vortex states are metastable, but also provides us with a 
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consistent QFT description of the unstable behavior. We examine some can- 
didates of states as to whether all the conditions are satisfied. In Sec. 8, for the 
physical states examined in Sec. 7, using Kubo's LRT, we calculate the density 
response of the system against the external perturbation of both impulsive and 
periodic types. As a result, we show that the complex eigenvalues correspond 
to the instability against some external time-dependent perturbations to the 
condensates. Section 9 is devoted to summary. 



2 Model Action and Hamiltonian 



We start with the following action to describe the trapped BEC of neutral 

atoms, 



d x 



i)\x){T -K- V{x) + fi}*P(x) - ^ 2 (x)i/j 2 (x 



where x — (x,t), and 



V(x) = ±m(u,lx*+u,y + uy), (4) 

with the mass of the neutral atoms m, the chemical potential fj, and the cou- 
pling constant g. Here we have written the trapping potential of harmonic 
type, but the discussions in this paper are valid for more general type of po- 
tentials. Throughout this paper h is set to be unity. 

In the terminology of the canonical operator formalism, let us divide the orig- 
inal field ift(x) into the classical and quantum parts as 

ij){x) = C(x) + <p(x) , (5) 

where it is assumed that the c-number function ((x), which is the order param- 
eter, is time-independent. Note that the function ((x) is essentially complex- 
valued if the condensate contains vortices, e.g., C{ x ) ~ e%Kd where 9 is an angle 
around a vortex line and n is a winding number. Equation (5) is substituted 
into Eq. (1), and the action is rewritten in terms of ((x) and (p(x) as follows: 

S = S + S 1 + S 2 + S 3A , (6) 

where 
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5 = J d'x [C(x){-K - V(x) + a*K(«) " f \C(x) 

5 1 = f*x{c{x) 



(7) 



-K-V(x) + ^-g\ax)\ 2 \ v{x) 

+ip \ x ) \- K - V{x) + /i — g\a*)\ 2 ] C(*)} , (8) 

5 2 = y rf 4 x jV(x) [T-K- V(x) + /i] ^(x) 

-f [4 |C(«)| V(x)¥»(x) + C*W (*) + C{x) V V{x)] } , (9) 

^3,4 = jd A x[-g [ax)^ 2 {x) V {x) + C{x)^{x)^{x)] -±<p*(x)<?(x)}. 

(10) 



At the tree level, the c-number function C( x ) satisfies 

K + V(x)-t, + g\((x)\ 2 ] C(x) = 0, 



(11) 



which is called GP equation [32]. The condensate particle number N c is given 
by 

N c = Jd 3 x\((x)\ 2 . (12) 

We move to the canonical formalism in the interaction representation. The 
total Hamiltonian of the system is now written as 



H = H + H mt + const. , 



(13) 



where 



H = jd?x\ <p\x){K + V(x) - fj}ip(x) 
9 



4|C(^)| 2 ^ t W^(^) + C\x)f{x) + C(x)^(x)} } , (14) 

H mt = Jd'x {g[c(x)$*(x)$(x) + C(x)<p\x)<?(x^ + ^ 2 (x)<p 2 (x)} . 

(15) 



It is emphasized that the CCRs, 

[<p(x, t), ] (x', t)\ = S 3 (x - x') 

and 

[<p(x,t),<p{x',t)\ = [fi(x,t),fi(x',t)} = 0, 
must hold in consistent QFT. 



(16) 
(17) 
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We introduce the new function f(x) which is defined as 



/(«) = ^(*) > ( 18 ) 

which implies that f(x) is normalized to be unity: 

Jd 3 x\f(x)f = l. (19) 

We can rewrite the GP equation (11) with f(x) as 

C(x)f(x)-M(x)f*(x) = 0, (20) 
where the following notations are introduced: 

C(x) = K + V(x) -fi + 2gN c \f(x)\ 2 , (21) 
M(x)=gN c f 2 (x). (22) 



3 Real eigenvalues of Bogoliubov-de Gennes Equations and Quasi- 
Particle Picture 

In this section, we review the BdG approach with real eigenvalues including 
zero-mode [33,34]. We consider the following set of equations for excitation 
modes, called the BdG equations: 

C{x)u n {x) - M(x)v n (x) = E n u n (x) , (23) 
C(x)v n (x) - M*{x)u n {x) = -E n v n (x) (24) 

for n — 1, 2, • • • . For the zero-mode, while f(x) satisfies Eq. (20), another 
function denoted by h(x) has to be introduced as 

C(x)h(x)+M(x)h*(x) = jf(x), (25) 

where J is a real constant which will be determined by the condition (28). Here 
we assumed that the phase part of h(x) is the same as that of f(x). This does 
not affect the generality of our discussion in this paper. The eigenf unctions 
u n (x), v n (x), fix) and h(x) satisfy the following orthogonal and completeness 
conditions: 
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J d 3 x [u* n {x)u n ,(x) - v* n (x)v nl {x)\ = 5 nn , , 
J d 3 x [u n (x)v n >(x) - v n (x)u n >(x)]=0, 

Jd 3 x [u* n (x)f(x)-v* n (x)f*(x)]=0, 
Jd 3 x [u* n (x)h(x) + v* n (x)h*(x)]=0, 



(26) 
(27) 
(28) 
(29) 
(30) 



and 



71=1 



£ [Un{x)u* n {x') - V* n (x)v n (x')\ + [f(x)h*(x') + h(x)r( X ')]=5 3 (x - X ') , 

(31) 

oo 

[u n (x)v* n (x') - v* n (x)u n (x')} - [f(x)h(x') - h(x)f(x>)} = . (32) 



n=l 



For real eigenvalues E n , it is known that the quasi-particle picture is obtained 
by diagonalization of the unperturbed Hamiltonian (14). The expansions of 
the quantum fields with respect to the annihilation- and creation-operators of 
the quasi-particle a n and d4 are given as 



= V(t)h(x) - iQ(t)f(x) + ]T [u n (x)a n (t) - v* n (x)ai (t)] , (33) 

n=l 

oo 

<p\x) = V(t)h*(x) + iQ(t)f*(x) + [<(x)6i(t) - v n (x)a n (t)} . (34) 

n=l 

The operators a n (t) and a4(£) are subject to the canonical commutation re- 
lations: [a n (t),al,(t)] = 8 nn i and [a n (t) , a n > (t)] = [a^t), a' n ,(t)] = 0, and Q(t) 
and V(t), which are called "quantum coordinates" (QCs), satisfy the canonical 
commutation relations: [Q(t),V(t)] = i and other commutations vanish. The 
relation between QCs and zero-energy particle mode is studied in Ref. [35]. 

Then, one can find the unperturbed Hamiltonian (14) is diagonalized in terms 
of a4(£), a n (t), Q(t) and V(t) as 



J oo 

H = —V 2 + Ena ] n a n + const. (35) 

ZI 71=1 

This determines the time dependence of the operators: 
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a n (t) =a n e~ 



iE n t 



«n(*) = & n e 



t JEnt 



Q(t) = Q + Vt, V(t)=V. 



(36) 



4 Doublet Notation 



We introduce the doublet notation, 



r(x) 



( ri(x) 
\r 2 {x) 



(37) 



Using this doublet notation, we define the following "inner product" for a pair 
of any doublets as 



(r,s) = J d z xr\x)a^s{x) 
= Jd 3 x (rt(x)r; 2 (x)) 



Sl{X) 

yO -1 J \s 2 (x) 
= J d 3 x [r\(x) s^x) - r* 2 {x) s 2 {x)\ , 

where (Tj represents the z-th Pauli matrix. 
We also define a (squared) "norm" of r as 



(38) 



| 2 = (r, r) = J d 3 x |ri(aj)| 2 — \r 2 (x 



(39) 



We note that this "norm" can be zero or negative due to the metric 03. 
The doublets of eigenfunctions of the BdG equations are defined as follows: 



x (x) 
x-x(x) 



\V n {x) ) 



(n = l,2,...), 



h(x) 
-h*(x) ) 



(40) 
(41) 
(42) 



We also introduce the symbol of 



y n (x) = aix* n (x) 



^ v*Jx] 



(n = l,2,-..) 



(43) 



Under the doublet notation, Eqs. (23) and (24) read as 

T(x)x n (x) = E n x n (x) , 

where 

T(x) = 



\ 



(44) 



C(x) -M(x) 
M*(x) -C{x) ) 
Using this representation we can rewrite the GP equation (20) as 

T(x)x (x) = 0, 

and Eq. (25) as 

T(x)x- 1 (x) = jx (x) . 

Note that ||a;o|| 2 — ll^-ill 2 — 0- The orthogonal conditions (26)-(30) and the 
completeness conditions are rewritten, respectively, as 



(45) 



(46) 



(47) 



(x n , X n ') $nn' i 
(y n ,yn>) =-8nn> , 

(y„,av)=0, 
(x ,x- 1 ) = l, 

(x n , Xq) = (x n , X_i) = , 

for n, n' — 1, 2, • • • , and 



(48) 
(49) 
(50) 
(51) 
(52) 



53 [ x n{x)xl(x') - y n (x) yi{x')\ + x (x)xl 1 (x') + x_ 1 (a;)xj(a3 / ) 

n=l 

= a 3 5 3 (x - x') . (53) 
The CCRs (16) and (17) are rewritten in the doublet expression as 

Mx, t),&(x', t)] = a 3 5 3 (x - x') , (54) 

where 



l>(x) 



' <f{x) N 
\$(x) J 



(55) 
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5 Complex Eigenvalues of Bogoliubov-de Gennes Equations 

Hereafter, we consider the situation in which the BdG equations have com- 
plex eigenvalues. Many numerical calculations [14,15,16,17,18] show that some 
complex eigenvalues appear in the presence of highly quantized vortices. We 
show some basic mathematical properties of the eigenfunctions, and then a 
symmetric property for the BdG equations. Finally we give a complete set 
including the eigenfunctions which belong to the complex eigenvalues. This 
plays a crucial role in constructing QFT involving the complex eigenvalues. 

5.1 Mathematical Preliminaries 

First, it is important to see that if E n is a complex number (lm(E n ) ^ 0), 
then, \\x n \\ 2 = 0. This fact is shown from the following relation: 

E n (x n , Xn) (Tx n , Xn) {x n , TXn) E n (x nj Xn) . (56) 

The same discussion of this fact is given in Ref. [15]. The similar analyses, 
related to this fact as well as the next statement shown below, are found for 
the special cases of the BECs with single multiply quantized vortex [16,18]. 

Next, let us see that if 

(x k , x rn ) ^ , (57) 

then, El = E m . Here x k and the eigenvectors whose eigenvalues are 

E k and E m , respectively. This is shown from 

E k \ x ki Xm) — (Txk, x m ) = (x k , Tx m ) = E m {x k , x n ) . (58) 

Note that x k and x m are not necessarily complex modes. 

We recall that the condensate can have complex eigenvalues when it has flow 
in a optical lattice or a highly quantized vortex. In those cases, the condition 
(57) holds [13,16,18], and the complex eigenvalues arise as a conjugate pair. 

As a contraposition of the above statement, we can see that if 

E* k ^E m , (59) 

then, (x k , x m ) = 0. This provides the sufficient condition for the orthogonality 
between different modes. In particular, two complex modes are orthogonal to 
each other if they satisfy the condition (59). We can also see that any complex 
mode is orthogonal to all real modes. 
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5.2 Symmetry for Bog oliubov-de Gennes Equations 



It is known that there exists a symmetry for the BdG equations when their 
eigenvalues are real [34] . For the BdG equations with complex eigenvalues, we 
find that if x n is an eigenvector which belongs to eigenvalue E n : 

T(x)x n (x) = E n x n (x) , (60) 

then, the following relation holds, 

T(x)y n (x) = -E* n y n (x) . (61) 

This implies that y n is also an eigenvector, belonging to the eigenvalue —E*. 
Note that the following relation also holds, 

\\y n \\ 2 = -\\x n \\ 2 . (62) 

To check the above relation, one uses c^TVi = — T. 

The similar analyses are found for the special cases of the BEC system with 
a single multiply quantized vortex [16,18] and of the BEC system under the 
periodic potential [13]. 



5.3 Complete Set with Complex Modes 



For QFT, we need a complete set in order to expand quantum fields. In this 
subsection, we present a complete set with complex modes. The Roman indices 
(n, m, ...) and the Greek ones (/x, v, ...), which are integers, are used for the 
real modes and for the complex ones, respectively, in our discussions below. 

Let us assume that there is a set of some pairs of complex modes, say {x^, x^}, 
which satisfy 



Tx M = E^ , (63) 
Tx^^ — E^^x*^ . (64) 

For different /x, it is also assumed that E^ is equal to neither E v nor E w for 
H 7^ v, meaning no degeneracy in complex modes. We adjust the normalization 
constants, so that the "inner product" is simply given by 

{%n i Xifv) = 5^ u . (65) 

The results of Subsec. 5.1 imply E 1 * = E^ and 

(v» , y*v) = -<V > ( 66 ) 



n 



where 



y*v = crix 



= <J\X 




(67) 
(68) 



For later convenience, we introduce the elements of and 



x 



x 




(69) 



(70) 



This way we can construct the complete set which includes the complex modes 
and is consistent with all the mathematical properties of the complex modes 
presented in Section 5.1 and Eqs. (65) and (66) as 



6 Representation of Free Hamiltonian and Expansion of Quantum 
Fields in Terms of Complex Mode Wave Functions 

6. 1 Representation of Free Hamiltonian and its Properties 

As we have the complete set with complex modes, we can expand the field 
operators. We discuss first in the Schrodinger picture and next move to the 
interaction picture. 

According to the completeness condition (71), the field operators are expanded 
in the Schrodinger picture as follows: 



[x n {x)x ] n {x') - y n (x)yl(x')\ + x (x)xl 1 (x') + rr_i(a:)4(0 



n 



+ J2 [ x ^ x ) x l^ x ') + x^(x)xl(x') - y li (x)yt li (x i ) - y^{x)y] i {x r ) 



= a 3 5 3 (x - x) . 



(71) 
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<f>( x ) = E [ & nU n (x) - & n v*(x)] + Vh(x) - iQ/(a>) 

n 

+ E [^W*) + B^{x) - Alv;(x) - B^vl^x) 



+ Vh*(x)+iQf*(x) 



n 



(72) 



(73) 



where the sets of operators {d: n ,a4} and are associated with the real 

modes and the zero-mode as discussed in Sec. 3, respectively, while the opera- 
tors Aft, At , 13^, and fit are newly introduced in connection with the complex 
modes. To be consistent with the CCRs, the new operators for complex modes 
have to satisfy the following commutation relations [15,22]: 



Ajuj 


si 




Afj,, 




= 0, 


Bf,, 


si 


= 0, 




B u 


= 0. 



(74) 
(75) 
(76) 
(77) 



Using the above representations, one can rewrite the free Hamiltonian H as 



(78) 



One can easily check the hermiticity of H in this representation. For later 
convenience, we denote the complex mode sector of H by H^: 



^ C) = E + E A\ 



(79) 



and we call Hq the complex mode sector Hamiltonian. One finds the time 
evolution of A^ and B^ as follows: 



A,{t) = e~ iE ^A 



4(t) 

B,(t) 
Bt(t) 



-iE*t 



B 



(80) 
(81) 
(82) 
(83) 



Note that this representation is consistent with the fact that H is time- 
independent. Thus the emergence of the complex eigenvalues does not imply 
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the instability of the system directly. This point will be discussed in detail in 
Sec. 7. 



Now, we obtain the quantum field in the interaction representation, 



<p(x) 



[a n e- tEnt u n (x) - ale^v^x)] + Vh(x) - i(Q + Vt)f(x) 

n 



-Ale^v;(x) - Ble^v^x) 



(84) 



£ [ale^uKx) - a n e- iE ^v n (x 

n 

+ E \A^%{x) + ^e^O) 

-A„e- iE »\(x) - B^e'^v^x) 
One can easily check that they satisfy the CCRs, 



+ Vh*(x) + i(Q + Vt)f*(x) 



(85) 



<p(x,t),<ft(x',t) =5 3 (x - x ) , 
[£(se, t),<p(x', t)} = W{x, t), fi(x', t) 



0. 



(86) 
(87) 



In order to study the complex mode operators and further, we represent 
them by the two kinds of bosonic operators [15]. We introduce the operators 
b/j,, frjj, and as follows: 



B,= ^(k-tbl), 
or equivalently, 



(89) 



The above definitions bring the following commutation relations: 



(90) 
(91) 
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= 8, 



(92) 
(93) 



with the other vanishing commutations. Similarly, one can easily check that 
the bosonic representations (88) and (89) satisfy the relations (74)-(77). 

Now, we can rewrite the Hamiltonian Hq C ^ in terms of and as follow [18]: 



(94) 



Let us see that is not diagonalizable using usual bosonic representation. 
Now we consider the Bogoliubov transformations, given by 



b\i — c n tn s ^\l ■> 
where the real numbers and s M satisfy 

2 2 1 

c M - s M = 1 , 

which ensures the bosonic commutation relations for ^ and 



(95) 
(96) 

(97) 



= 5 



(98) 
(99) 



and zero for the other commutations. Then is rewritten in terms of £ and 
I as 



= J2 {ME,) - 2c^lm(E,)}il 4 - {Re(^) + 2c^Im(£ M )} |J £„ 



+ Im(^)(<£ + + Im(^)(cJ + ^)U, 



(100) 



The Hamiltonian could be diagonal only if the following relation for each 
H were fulfilled, 

lm(E^(cl + sl)=0. (101) 

But from Eq. (97), we have cj + sj = 2sJ + 1 > 1, and Eq. (101) would 
lead to In^E^) = for each fi. Thus when the eigenvalue is complex, the 
Bogoliubov transformation can not diagonalize Hq\ We therefore conclude 
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that the complex mode sector Hamiltonian Hq is not diagonalizable in the 
bosonic representation. 

* (c) 

6.2 Eigenstates of Complex Mode Sector Hamiltonian Hq 



In this subsection, we construct eigenstates of Hq C \ starting from the vacuum 
of 6 M and 6 M . 

First, we consider a transformation from the set of the bosonic operators 
{6^, 6^, 6 M , 6^} to the set of the complex mode operators {A^, A^, B^, Bj^}. 
One finds the following explicit transformations: 



IX IX IX 



and 



where 



= exp 

Note that the operator is not unitary but hermitian, 

/x vy A* ' 

and its inverse operator is given as 



W~ l = exp 



.7T 



For later convenience, we introduce the notation of 



w- 1 = ® w- 1 . 



(102) 
(103) 



(104) 
(105) 

(106) 
(107) 

(108) 



(109) 
(110) 



Next, we construct states which are annihilated by A^, A^, B^ and B^. In 
order to find those states, we introduce a state |0)) defined by 

|0» = ®|0»„, (HI) 
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where |0)) is the tensor product of the vacuum of 6 M and b^. 

|0)) M =|0) b ,®|0)^, 

where 

Obviously, the state |0)) is annihilated by both ft M and 6 M , i.e., 

/ ^|0» M = 0, 

Now, we introduce the states |0) A and |0) B defined by 
|0) A = W|0», 

|o> fl =w- 1 |o», 

Because of hermiticity of (107), we obtain ^(0| and B (0| as 

a(0\ =((0\W, 
B (0\=((0\W- 1 . 

One can easily find the following relations by using Eqs. (102)-(105) 

A II \0) A =0, 
^|0>a= 0, 
B^\0) B =0, 
BI\0)b = 0, 

and 

a<0|4=0, 

b<o|4 = o. 

B<0|St = 0, 

17 



which are valid for arbitrary fx. Hereafter we call \0) A and |0) B zero states of 
and B^, respectively. 

These zero states are eigenstates of Hq C \ 



Ht ] \0) A =- 



^ c) |o), 



UK 



|0)a, 
\0)b- 



Note the following relations for the bra states: 



(127) 
(128) 



a(0|^ c) 

*<o|^ c) 



T,e; 



a(0\, 
b(0\. 



(129) 
(130) 



Using Eqs. (115) and (116), one easily derives 

a(0\0)b = b(0\0)a = 1- 



(131) 



It is important to notice that both of a(0|0)^ and b(0|0) b diverge. To show 
this fact, we rewrite ^(010)^ in the bosonic representation: 



A<o|o> A = n[4o|^ 2 |o>> 

=nL((oie 

Define a new state \9) by 
\0) = 



^)|0», 



(132) 



exp(- In cosh 9) exp tanh 6>) (0)), 



(133) 



where the constant term on the right-hand side is determined by the normal- 
izatoin condition (^l^) = 1. One finds the limiting behavior 



«0|<9> -> oo (e^-i 



(134) 
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which implies that a(0|0)a is divergent. Similarly, it can be proven that s(0|0) b 
is also divergent. 



6. 3 Complete Set of States for Complex Modes 



We start from 



with the definition 



iB E 



(135) 



(136) 



( c ) 

Here, we defined the new symbol = which is an equal sign only for parts 
belonging to the complex eigenvalues, e.g., H = ■ 

Having W and W' 1 act on Eq. (135) from the left and right, respectively, we 
obtain 

1 ( = } \0)a b(0\ 

+ Y,H\ Q )ab(0\A,]-Y,[b»\0)ab(0\AI] 



+ 



(137) 



which is the complete set of states for the complex modes. Performing the 
similar procedure, one obtains another representation: 

1 = \0) B a(0\ 



+ 



£ [Al\0) B A (0\B»] - £ [A,\0) B A (0\Bl 



1 K 2 \*)b a(0\BI 



+ 



(138) 



Important observations here are that (i) the "natural" conjugate of \0) A is 
b(0\ and vice versa and that (ii) the minus sign enters the complete set. The 
latter originates from the fact that the state space is an indefinite metric one. 



19 



7 Physical States 



In the previous section, we studied the properties of H when there exist the 
complex eigenvalues. It turned out that the state space is not the simple Fock 
one but the indefinite metric one. Therefore, we need to impose appropriate 
conditions to construct a restricted physical state space. The physical states 
are required to reproduce the current experimental results. It was revealed in 
the BEC experiment [8] that the condensate with a doubly quantized vortex is 
metastable and that its life time is not related with the value of the imaginary 
part of the complex eigenvalue of the BdG equations directly. Note that, in 
QFT, even if the system is unstable, we start with a static picture and then the 
unstable behavior is described. Recall for example the description of the Beli- 
aev process [36], that is, we first construct the static representation with the 
stable unperturbed Hamiltonian and then the decay processes are described as 
the higher order of perturbation. Here we take the density distribution as the 
static quantity in accordance with the experiment [8]. As discussed later, it is 
expected that the unstable behavior occurs due to the external perturbation 
in this case. 

We now require the following four physical state conditions (PSCs): 

(i) (*\j,(x)\*) = {(x), 

(ii) {^l^ (x)ij}(x) |\&) is time-independent, 

(iii) (\&|G f |\I>) is real, when G is an hermitian operator, 

(iv) (tf |tf> = 1, 

where |\&) is a physical state and (^\ is its "natural" conjugate (see the com- 
ment below Eq. (138)). We call |\&) and (^\ physical states if they satisfy the 
four conditions above. The first and second conditions mean that the order 
parameter and the density distribution are static. The third condition ensures 
that the expectation value of any hermitian operator can be interpreted as 
physical quantity. The fourth condition is necessary for the probabilistic in- 
terpretation. Hereafter, we assume for simplicity that there is only one pair of 
complex modes, and denote their eigenvalues, eigenfunctions and operators by 
E and E*, u(x), v(x), u*(x), and i>*(cc), and A, A*, B, and W , respectively. 

From the equation ip(x) = ((x) + (p(x), the first condition of the PSCs (i) is 
equivalent to 

(i') (tf|0(x)|¥) = O. 

If and (^1 satisfy the condition (i'), the second condition (ii) states that 
(ii') (\I'|(^" l "(a;)( l 3(x)|\E') is time- independent. 
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From the discussion above, one can confirm whether a state is physical one or 
not, by checking first the condition (i'), next the condition (ii') and (iii) , and 
finally the condition (iv). Let us examine some states, following these steps. 

7. 1 Candidate 1 — Vacuum of b and b 

We consider the vacuum |0)) as the complex part of the first candidate for the 
physical state namely, 

l*>i = |0»- (139) 
The complex part of the corresponding conjugate state is ((0|, i.e., 

= «0|. (140) 
One can easily see from the following relation: 

«0|i|0» = «0|it| » = « |B|0» = «0|it| » = , (141) 
this state satisfies the first condition (i'). 
We move to the condition (ii'). It is easy to obtain 

«0|iit| 0>> = (( |SSt|0» = \ . (142) 

Combining these results with 

<p\x)(p{x) 

i iit e 2(ims)t [ u ( aJ ) u *( a .) + v ( x ) v *( x )] 

- AAe~ 2[Et u(x)v(x) - A^A^e^u^xyix) 
+ B&e- 2 ^ 1 [u,(x)v*(x) + v*(x)vl(x)] 

- BBe- 2iEH u*{x)v*{x) - &&e 2iEt ul{x)vt{x) 

- ABe- 2iiKcE)t [u(x)v*(x) + u*(x)v(x)] 

- A^&e 2i ( ReE » [u*(xX(x) + ut(x)v*(x)} 
+ A ] Bu* '(x)u*(x) + BA ] v*(x)v*(x) 

+ &Au(x)ul(x) + AB*v(x)vl(x) , (143) 

we can see that ((0\tp' (x)tp(x)\0)) is time dependent. This fact means that 
the condition (ii') is not satisfied. So the pair of the states and is 
excluded from physical states. 
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7.2 Candidate 2 — Zero States 



In this subsection, we examine the zero states (see the definition below Eq. (126)) 
as the second candidate for the physical states. 

Let us pick up a pair of \0)a and its "natural" conjugate b(0|, namely, we 
choose the complex part of the second candidate for the physical states |^) 2 
and 2 (^| as 

|*> 2 = |0)a, 2(^1 = b(0\ (144) 
(see again the comment below Eq. (138)). We can easily check the following 
relations, 

B (0\A\0) A = B (0\A^\0) A = 0, (145) 
B<0|B|0) A = B (0|B t |0) A = 0, (146) 

which means that the condition (i') is satisfied. 
Next, using Eq. (143), we obtain 

2(*\<p\x)<p(x)\y) 2 ( = } -u*(x)u*(x) + v{x)vl(x). (147) 

This shows that the condition (ii') is satisfied but the condition (iii) is not. 

Thus we conclude that this pair of the zero states is excluded from physical 
states. Similarly, the other pair of |0) s and ^(0| doesn't satisfy the PSCs 
either. 



7.3 Candidate 3 — Direct Sum of Zero States 



In this subsection, we consider the direct sum 

\+)^^m A ®\0) B ) (148) 

as the complex part of the third candidate for the physical state |^) 3 , namely, 

l*>3=|+>. (149) 
We may introduce such a direct sum because these two zero states belong to 
the different complete sets, as one can see from Eqs. (137) and (138). Note 
also that in this case the complex part of the corresponding conjugate state 
is given as 

*m = {+\ = ^(B(0\®A{0\). (150) 



22 



Using Eqs. (145) and (146) and their complex conjugates, one can see that 
the condition (i') is satisfied. 

Next, using Eq. (143), we obtain 

3 m^(x)0(x)\^) 3 ( = } -Re (u*(x)u*(x) - v(x)v*(x)) . (151) 

This time, the conditions (ii') and (iii) are both satisfied. In particular, note the 
fact that the number density is time-independent. Thus, as already mentioned 
in Sec. 6, the emergence of the complex eigenmodes does not immediately 
imply the instability of the system in our approach. It is easy to show that the 
condition (iv) is also satisfied. So, we can conclude that the third candidate 
for the physical states \^/) 3 and 3 (\l/| satisfy all the PSCs. 



8 Linear Response for Direct Sum of Zero States 

So far, we have developed the description of QFT in the static picture when 
the BdG equations provide complex eigenvalues. Then the complex eigenvalues 
are not directly connected with the instability of the system. If we choose the 
direct sum of the zero states as the physical states, the distribution of the 
quasi-particle is stable. It is not yet clear how the complex modes give rise 
to any instability. In order to clarify this point, we consider the response 
of the condensate against the time-dependent external perturbations, whose 
expression is derived from the linear response theory [28,29,30]. It will turn out 
that the response is remarkably affected by the presence of complex modes. 

We consider the following external perturbation which corresponds to time- 
dependent modification of the trap: 

H cx (t)= Jd 3 x^(x)5V cx (x)^(x) . (152) 

From the linear response theory, the response of the condensate density dis- 
tribution 5(p(x,t)) is given as 

5(p(x)} = Jd 4 x'G R (x,x',t-t')SV cx (x'), (153) 

where 

p(x) = ft(x)4>(x) (154) 
and 5V e x(x) is a time- dependent external potential. The retarded Green func- 
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tion Gr(x, x', t — t') is expressed as 

G R (x, x', t-t') = -iO(t - t)([ft(x)4>(x), ^(x')ij(x')]) 
= -iC(x)<(x')9(t - (0(x, t)<p\x', f )) 

+ iC(x)c(x')e(t-t')(^(x',t') ( p(x,t)) 
- i((x)C(x')e(t - 1') ($(x, t)<p(x', 0) 

+ <C(aj)r(*')«(*-* , )^(« , ,< , )^ t (a',*)> 

+ -.., (155) 

where the symbol "• • •" stands for the higher order terms of ip and </?', and 
(•) represents the expectation value of some physical states. In the following 
discussion, we consider the case in which there is a single pair of complex 
modes as in the previous section, and employ the states |\&) 3 and 3 (\l/| as the 

physical states, whose complex parts are given as |\&) 3 = |+) and 3 (^| = (+| 
(see Eqs. (148)-(150)). 

The type of the time-dependent external potential should correspond to an ex- 
perimental setup, and we choose two types: One is an impulsive force, and the 
other is an oscillating one. In the case of an oscillating external perturbation, 
we assume that the perturbation starts at to and is kept through. 

First, we focus on the impulsive external potential. In this case, the external 
potential is given as 

SV ex (x) = S(t - t )SV(x) . (156) 
From Eq. (155), one obtains the following density response, 

5{p( x )) = -iC(x)u(x) (^J<Px' C(x')uZ(x')5V(x')) e - lE ^ to) 
-i(*(x)u*(x) (Jd 3 x'((x')u*(x')5V(x'^j e-^-' ) 
+ iC(x)v*(x) d 3 x' C(x')v,(x')5V(x')j e <B '(*-*o) 
+ iC(x)vl(x) (Jd 3 x'((x')v(x')5V(x')^ e iE ^ 
- i((x)v(x) (I d 3 x' C(x'K(x')5V(x')j e' tE{t ~ to) 
-i((x)v*(x) (Jd 3 x'C(x')v*(x')SV(x'^j e-^-' ) 
+ i((x)u*(x) (| d 3 x' C(x')u,(x')5V(x')j e tE ' t{t ~ to) 
+ K{x)ul{x) (| d 3 x' C{x')u{x')5V{x')} e iEit ~ tQ) . (157) 

It is seen from the above expression that the density response S(p(x)) blows 
up or damps in general if the complex modes exist. In this sense, we can 
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say that the complex eigenmodes cause the instability against the external 
perturbation. This result is consistent with that of the analyses using the 
TDGP equation [17,18]. 



Next, apply an oscillating external potential which is given as 

8V ex (x) = 9(t — t )R(x) coscut . 



(158) 



Here R(x) is some function of x representing a type of oscillation. For example, 
one may choose R(x) = x in the case of the dipole oscillation of the condensate 
to the x-direction. The general expression of the response of the condensate 
is given as 



5(p(x)) = - 



ik[{, u ; 



■C(x)u(x)e lMo ( e ^*-*°) - e - lE{t - to) ) 



2 uj + E 



-iE*(t-t ) 



) 



2 oj-E 

- - K[C,U * ] C(x)u*(xy Mo f e -(*-*o) 
2 Cu + E* s v ; v ; V 

1 K\^v4 iMo , Mt _ to) _ lE * {t . to) \ 

2 Cj - E* y ' v ; V ) 

- - K ^ ,V *k *(x)v*(x)e-^ t0 f e --(*-*°) - e 1 ^^) 
2 u + E* y 1 y 1 V J 

~ l^4c(x)v:(x)e ii:,t0 (e^-V - e iE «- t0 A 
2 cj — E ^ ' 

ik[{,v 



+ 



2 Co + E 
2 uj + E 



+ 
+ 



2 u-E 
l K[C,v' 
2 uj + E* 
I K[C,V 
2 uj-E* 
l K[C,u, 
2 lj-E* 
lK[C,u, 



C(xK(x)e- lM " (e-^-'°) - e iS(t - t()) ) 

^-C(x)v(x)e-^ t " (e-^ 1 -^ - e - iE{t - to) ) 
k(x)v*(x)e lM » (e iCj{t - to) - e - tE * {t - ta) ) 
k(x)v4x)e- iM " ( e -^(*-*o) - e'^* (t - to) ) 
lc(a;)^(a;)e^ ( e ^-*°) - e ^* ( '~' o) ) 



2 £ + £ 



(159) 



25 



where 




(161) 



(160) 



This again shows that the instability against the external perturbation comes 
from the complex eigenmodes. 



9 Summary 

In this paper, we have constructed the description of QFT (quantum field the- 
ory) starting with a static picture and have formulated Kubo's LRT (linear 
response theory) in the case where complex eigenmodes of the BdG equations 
arise. First, we have summarized mathematical properties of the eigenfunc- 
tions. These properties guide us to a consistent form of the complete set. Using 
the complete set, we expand quantum fields and obtain the representation of 
the free Hamiltonian. We have shown that the Hamiltonian is not diagonal- 
izable in the conventional bosonic representation. The eigenstates of the free 
Hamiltonian can be constructed, although their "norms" diverge. These in- 
tricate circumstances, when the complex eigenmodes are involved, come from 
the fact that an indefinite metric space comes in. To deal with an indefinite 
metric and to construct a physical state space, we have introduced the PSCs 
(physical state conditions), reflecting the metastability of doubly quantized 
vortex states of BEC. We investigate the candidates for the physical states. 
It turns out that if we employ the vacuum of b and b (Candidate 1) and the 
zero states \0) A or |0) B (Candidate 2) as the complex part of the states, they 
don't satisfy all the PSCs. But, if we choose the direct sum of the zero states 
(Candidate 3) as the complex part of the states, they satisfy the PSCs. Using 
the direct sum of the zero states as the complex part of the physical states, we 
can start with the stable representation in QFT and calculate the (retarded) 
Green function. Finally the density response against the external perturba- 
tions (the impulse and oscillating types) is derived. This result shows blow-up 
and damping behavior of the fluctuations, which is qualitatively consistent 
with the result of the analyses using the TDGP equations [17,18]. Thus the 
complex eigenmodes cause the instability against the external perturbation. 

We mention that the higher order terms of ip and tp' have contributions to 
the density response, and have the poles at E n ± E and so on, which are the 
contribution beyond the calculations based on the TDGP equation. 
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